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Abstract 

In this article we study a simple spin model which has a non-interacting Hamiltonian but 
constrained dynamics. The model, which is a simplification of a purely toplogical cellular 
model displays glassy behaviour, involves activated processes and exhibits two-step re- 

laxation. This is a consequence of the existence of annihilation-diffusion processes on two 
distinct time-scales, one temperature independent and the other an exponential function of 
inverse temperature. In fact, there are several such inter-coupled microscopic processes and 
great richness therein. Two versions of the model are considered, one with a single absorbing 
ground state and the other with a highly degenerate ground state. These display qualita- 
tively similar but quantitatively distinct macroscopic behaviour and related, but different, 
microscopic behaviour. 



1 Introduction 

In a recent paper Davison and Sherrington focussed on a purely topological tiling model Q 
which exhibited glassy dynamical behaviour. This was driven by a desire to investigate the be- 
haviour of supercooled liquids in a model which contained as few parameters as possible, but still 
displayed the relevant physics. In this paper we continue the minimalistic approach, and consider 
systems which are conceptually similar to the topological model, but have the major advantage of 
involving variables (spins) based on a fixed lattice. 

The original model was that of a two-dimensional topological froth, constructed by tiling the 
plane using three-fold vertices only. However, unlike soap froths, the energy was determined by 
the deviation of the cell topologies from a perfect hexagonal tiling: E = Yli i n i ~ 6) 2 , where 
rij is the number of sides of cell i, and the system evolved solely through stochastic Tl-micro- 
dynamics with Glauber-Kawasaki probabilities. We found that there was two-step relaxation 
at low temperatures, and formed the following conceptual picture describing the evolution of 
the system. At low temperatures the system consists mainly of six-sided cells, and there are 
two processes dominating the behaviour: on a fast time-scale, pairs of pentagon-heptagon defects 
diffuse freely through the hexagonal background, and on a slower time-scale isolated defects absorb 
or create pairs of defects. The latter can be considered to be an activated process, as it is 
energetically unfavourable for such a pair of defects to be created. This conceptual framework 
of two different processes leading to both fast and slow dynamics is directly applicable to the 
lattice-based spin model with which this paper is concerned, and thus one might expect it to yield 
qualitatively similar results to the topological model. Given that this spin model is computationally 
simpler and more tractable, one might hope to be able to probe more deeply, and to investigate 
to what extent one finds the same features displayed by other kinetically constrained models 
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(a) The Tl moves for the topological froth. 




(b) The spin-flip rules for the present model: the total spin is conserved and the choice of whether 
to attempt the upper or lower signs is made randomly at each time-step. The quadruplet of cells 
is identified uniquely by the dashed edge e. 

Figure 1: Spin- flip rules for the 3-state spin model and the topological froth. 



(for examples see p|-|Tl]] ) . In fact, as we show below, these expectations are borne out and it is 
possible to provide a physical understanding of the features observed. It is also both possible 
and instructive to consider a generalization with a very different (highly degenerate) ground state, 
which we cover in the latter half of the paper. 

2 The Model 

The model we use comprises a perfect hexagonal tiling of the plane, with a variable (spin) associ- 
ated with each hexagon. The spin in cell i, denoted by Si, is restricted to the values 0, ±1. This 
is analogous to a topological froth in which the cells are all pentagons, hexagons or heptagons, so 
that the topological charge qi = (6 — of each cell is restricted to the values 0,±1. However, 
the topological froth model has a non-periodic and dynamically changing cell structure, whereas 
this model is firmly fixed on a perfect hexagonal tiling. 
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We define the energy as follows: 



N 

E = D^ (1) 

i=i 

where N is the total number of cells/spins in the system, and D may be positive or negative. 
The case of D > emulates the original topological model where hexagonal order is energetically 
preferable, but D < is also of interest, as we shall show. The dynamics conserve the total spin 
of the system and we choose our starting configuration such that the total spin is always zero i.e. 

JY 

E s * = ° ( 2 ) 

i=l 

in analogy with the Euler rule {6 — rii) =0 which applies to a froth j fl2| , |l3] ]. Therefore the 
ground state for D > consists of Si — for all spins, whereas for D < the ground state is 
degenerate, with half the spins taking the value +1 and half taking the value —1. 

The system is allowed to evolve through local spin-flips which are similar to the Tl moves 
performed on topological froths. In the latter two adjacent cells have their topological charges 
decreased by 1, and their two common neighbours have their topological charges increased by 1: see 



Figure 1(a) . In the present model the allowed move-sets consist of choosing a pair of neighbouring 
cells and either increasing their spins by 1 unit, and decreasing those of their common neighbours 
by 1 unit, or vice versa; both possibilities need to be allowed to avoid chirality inhibiting movement 
of spins throughout the system. 

The actual dynamical process is as follows. At each time step an edge is chosen randomly on 



the hexagonal lattice; this defines a set of four cells as shown in Figure 1(b). A choice is then 
made randomly of whether to consider increasing the spins of the adjacent cells y and w (and 
thus decreasing v, x) or vice versa, with equal probabilities for both cases. The probability of 
actually performing the move is dependent on the energy change that would be incurred, and is 
given by a temperature-dependent Metropolis-Kawasaki^] algorithm. Specifically, assuming w and 
y have been chosen as candidates for an increase in spin, the energy change associated with these 
spin- flips on spins s v ,s w , s x , s y is: 

AE(s w , s v \ s v , s x ) = 2D(2 + s w + s y - s v - s x ) (3) 

and the probability P of actually performing this move is: 

P(s Wl s y ; s v ,s x ) = (1 - S Swt i)(l - S Sy> i)(l - 5 Sv -i)(l - 6 Sl -i) Mm[l, exp(-fiAE(s w , s y \ s v , s x ))] 

(4) 

where (3 is the inverse temperature. The ^-functions ensure that the spins are forbidden to take 
values other than ±1 or 0. 

The simple form of equation ([!]), with no interaction between the cells, shows that this system 
is thermodynamically trivial in equilibrium and all the static equilibrium properties are readily 
calculable. However, the microscopic dynamics are constrained and non-trivial, involving several 
spins simultaneously; this leads to glassy macro-dynamics. 

Most of the data results from simulations on a system of size N=9900, although, in order to 
perform more accurate fits, in certain cases the system size was increased to 160000. However, 
unless otherwise stated, one should assume the former system size is in use. Periodic boundary 
conditions are enforced in all cases. 

This paper is structured in the following manner: first we present results for D > 0, including 
a discussion on the processes involved in relaxation of the system; this is followed with results for 
D < 0. Finally there is a more general discussion. 

1 We choose to use the Metropolis algorithm rather than Glauber dynamics as in previous work because the 
qualitative features of the results show no dependence on which of these algorithms we choose, and Metropolis is 
the faster of these two. The reference to Kawasaki is included to emphasise the fact that although more than one 
spin is flipped at once, the total spin is conserved. 
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Brief review of the topological froth model 

As one of the aims of this paper is to show that the behaviour of this model is indeed qualitatively 
very similar to that of the topological froth, we shall briefly review the results for the topological 
froth before proceeding 

In simulations in which the system is cooled at a variety of different rates, one finds strong 
dependence of the energy on the cooling rate, with the system unable to attain equilibrium within 
any reasonable time-scale at very low temperatures. Measurements of a two-time auto-correlation 
function show evidence of two-step relaxation, with plateaux developing as the temperature is 
reduced. There are also clear signs of aging when the system is not in equilibrium. When the 
equilibrium correlation functions are rescaled by a suitably defined relaxation time, one finds that 
they collapse onto a master curve in the late-/3 relaxation regime, and that this master curve can be 
fitted by a von Schweidler law as predicted by Mode-coupling Theory (MCT) [Q. The relaxation 
time is well-described by an offset Arrhenius law, indicating strong glassy behaviour |L5| . Although 
MCT in fact predicts a power law, this does not fit the data particularly well, and neither does 
the oft-used Vogel-Fulcher law. 

By measuring a suitable temporal response function and plotting it parametrically against the 
appropriate correlation function (starting from a non-equilibrium configuration), the Fluctuation- 
Dissipation Theorem (FDT) is found to be upheld for times somewhat longer than required for 
the onset of the correlation function plateaux. For longer times FDT ceases to hold. After it is 
broken, one observes non-monotonic behaviour, a feature which has also been noted in several 
other models which can be considered to involve activation over energy barriers 

3 D > 

3.1 Relaxation Dynamics 

It is a simple matter to calculate the equilibrium behaviour of J| , for which one finds (in units of 

D = iy. 

E 2exp(-/3) 

iV l + 2exp(-/3) ' [ ' 

One has instant access to equilibrium states, as one can randomly place the appropriate number of 
±l's throughout the system to access a particular temperature. Using this feature, we may study 
the behaviour of the energy with temperature at a number of different cooling rates, starting from 
equilibrium at any chosen temperature. Above T = 1 the system equilibrates very rapidly at even 
the fastest cooling rate, so we have chosen a starting configuration of T = 1, and then cooling is 
carried out by waiting a time t w = jN at each temperature decrement of ST — 0.05. The results 
(averaged over 3 runs), with the equilibrium curve, are shown in Figure 2(a)| ; the system exhibits 



strong dependence of the energy upon the cooling rate. This is characteristic of glassy systems, 
and qualitatively similar to the results found for the purely topological froth. 

It is instructive to study also the behaviour of the system when subjected to a rapid quench 
from an infinite temperature {(3 = 0) configuration to a temperature at which it is allowed to 
evolve for a time t w . Figure |2(b)| shows the results from such a quench for a range of different 
temperatures. For longer values of t w we see a minimum develop - this strongly suggests activation 
is present: at very low temperatures the system cannot overcome the energy barriers and thus 
cannot access lower energy states. The temperature at which this minimum occurs is dependent 
on the waiting time t w . At very low temperatures, even at the largest waiting times employed 
(t w = 1000007V) the system is unable to reach energies below f - 0.09. 

One sees the significance of this value if we turn to the temporal behaviour of the energy. 



Figure |3(a)| is a plot of M against time, quenched from a fully disordered starting configuration 



corresponding to (3 — to the temperature in question. The initial decay of the energy is fast and 

E_ 

N 



independent of temperature until jt ~ 0.09, at which point one sees the existence of a plateau; 
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(a) The behaviour of the energy with temperature under slow cooling (D > 0). 
The values of t w are the waiting times at each point. 



T 




Temperature 



(b) The behaviour of the energy with temperature after a rapid quench (D > 0). 
The values of t w are the times, subsequent to the quench, at which the energy is 
measured. 

Figure 2: Energy against temperature for slow cooling and rapid quench. 
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(a) Energy against time. The curves for /3 < 5 reach their equilibrium values, whereas 
the others do not. 
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(b) Energy against Tint (where t is measured in units of N). 

Figure 3: The behaviour of the energy with time. 



(a) Annihilation of two dimers. This is also possible (b) Annihilation of a dimer through interaction 

if each spin is multiplied by -1. with a defect. The defect is shifted in position. 




(c) Free diffusion of a dimer in a background of zero 
spins. 



Figure 4: The dominant moves through which this model evolves. 



the time spent on this plateau is clearly dependent on temperature. Upon departure from the 
plateau, the energy relaxes directly to the appropriate equilibrium value; on the graph shown, 
the curves for [3 < 5 equilibrate within the time-scale of the simulation, whereas those for (3 > 5 
do not. This two-time behaviour, with one time-scale temperature independent and the other 
increasing with inverse temperature, can be understood as follows. The mechanism for lowering 
the energy is the annihilation of pairs of adjacent +1 and —1 spins (we shall refer to a +1, —1 pair 
as a 'dimer'). This occurs in two ways: (i) two n eighb ouring conjugate dimers can destroy each 



other to leave four zero spins as shown in Figure 4(a) , with a reduction in energy of 4 units, or 



(ii) a ±1 can annihilate a dimer, thus shifting its position and reducing the energy by 2 units as 
in Figure 4(b) ; note that the ±1 can be part of a dimer at a different orientation. Any such local 
arrangements present in the starting configuration will be eliminated quickly without need of any 
thermal excitation. Furthermore, dimers can move freely through a background of zero spins as 
shown in Figure 4(c) until they reach a local environment which favours annihilation, such as those 
previously mentioned. This diffusion occurs on a time-scale of 2 steps per spin, as one can alter 
the configuration on the left of the arrow in Figure 4(c) in 2 ways, one of which will be possible 



and one of which will be forbidden through the 5- functions in equation (Q). Consequently, the 
initial fast decay of the energy is temperature-independent, and of a diffusive character with an 
underlying time-scale of 2. This fast, diffusive process describes the behaviour of the energy until 
the plateau is reached. To describe it further, one has to consider isolated defects i.e. isolated 
spins of ±1. To remove these defects, they must be paired up with a conjugate isolated defect to 
form a dimer which can then diffuse freely as in Figure 4(c) , and eventually annihilate as in Figure 
4(a) and 4(b). An isolated defect can move through fortuitous collisions with existing dimers, but 
after the initial fast decay these dimers beco me rar e. Alternatively, a defect can move by creating 
a new dimer (reversing the arrow in Figure 4(b) ), at an activation energy cost of 2 units and 
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with a probability that scales as e -2 ^. As is clear from this figure, one can interpret the resulting 
configuration as a dimer-plus-defect in two ways; either of these two possible dimers can diffuse 
away freely if adjacent to two zero spins, or annihilate if adjacent to another dimer. This factor 
of 2 will cancel that introduced by the diffusive time-scale of the dimers. Thus again one has a 
diffusive process leading to a final reduction in energy, in this case with a time-scale of e 2/3 (since 
the time-scale for annihilation of the dimer is negligible compared to that to produce the dimer 
for e 2/3 3> 1)- Thus E(t)/N is expected to consist of two diffusive processes: a fast process of 
time-scale 2 decaying to a state of isolated defects, and a slow process of time-scale e 2/3 decaying 
to the equilibrium configuration. 

The plateau in E{t)/N can be seen more clearly if the time axis is rescaled to Tin t as in Figure 
[3(b)| ; as the temperature is decreased (/? increased) the curves tend to a sharp staircase form. This 
is reminiscent of the results found under such rescaling for other kinetically constrained models 
H [|; in these other models one observes several plateaux corresponding to several characteristic 
activation energies, but in this particular case the situation is simpler as there is only one dominant 
characteristic activation energy. 

Both the fast dimer-dimer annihilation and the slow defect-antidefect pairing are of the type 
usually designated as A + B — » 17 In the fast process, A and B are dimers and 'anti-dimers' 
i.e. a (+1,-1) dimer annihilating with a (-1,+1) anti-dimer; for the slow process, A and B are 
isolated defects of opposite sign. In detail these diffusion processes are more complicated than 
simple diffusion, but for A + B — > processes the standard asymptotic behaviour of the density 
is of the form t~ 3 , where d is the dimensionality. Therefore we suggest the same asymptotic 
(t/r)~ - 5 behaviour for each process^] and fit the following form to the energy: 

2 In fact, whilst the slow processes are isotropic, this is not the case for the fast dimer diffusion: the latter involves 
zig-zags at 30° to the axis perpendicular to the common edge between the two cells constituting the dimer. 
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Figure 6: The behaviour of E/N after a T = quench from equilibrium. 



E(t) [2 \ ( t\ , J t , 

- ' - W ' ] ' (a-e eq ) +e eq (6) 



where a is the plateau value, e eq is the energy per spin in equilibrium, and we expect both b 
and c to be approximately 0.5. The results are shown in Figure ||; it is clear that these fits are 
extremely good. In principle, the plateau value a can be calculated but here we note only that, 
as required, a is less than 0.25, which is the maximum energy for a T — frozen state. It is also 
less than the value corresponding to randomly removing dimers from the initial configuration to 
leave only isolated defects, which gives E{t)/N ~ 0.2; in fact, one finds a ~ 0.085 due to the effect 
of singleton pairing by dimer collisions. Further investigation of the characteristic time-scales is 
given below, in connection with the correlation function. 

The reader will note that the fast dimer absorption by an isolated defect has been neglected 
in the above fit; this can be characterised as type A + C — > + C, and standard asymptotic 
behaviour of the density for such a process is that of a stretched exponential ]l8|-pl[. There is 
also a move-set which involves two dimers interacting 'off-centre', such that one of the dimers and 
only one of the defects in the other dimer are altered. This leaves either two isolated defects or a 
dimer. However, in view of the excellent quality of the fit we do not consider these explicitly at 
this stage. 

We have also investigated the behaviour of the energy if one quenches to T — from an 
initial equilibrium configuration corresponding to a finite temperature Tj. This is in some sense 
an investigation of the inherent states of the system |p2[-pq|. At very low temperatures, the 
equilibrium state for Tj consists mainly of isolated defects - thus when quenched to T = the 
system very quickly reaches the inherent state. However, at higher initial temperatures this takes 
an extremely long time to happen. Thus we show in Figure || the energy measured at waiting 
time t w after a quench from the equilibrium configuration to T — 0. At low Tj the energy stays 
on the equilibrium curve, as there are no energetically favourable moves to be made. However, 
for higher Tj the energy tends towards a constant value given by the plateau in Figure ||. Thus 
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we see that if one quenches from an equilibrium configuration with energy above the plateau, 
the lowest energy one can reach is that of the plateau - to decrease energy further would require 
activated processes, which one cannot perform at T = 0. If one quenches from an equilibrium 
start-point with energy below that of the plateau, one cannot decrease the energy by much as 
there are very few energetically favourable moves to be made, and so the curves in Figure |6| do not 
deviate far from the equilibrium curve. The temperature at which one sees a crossover between 
these two types of behaviour is T ~ 0.35: this is the temperature at which activated processes 
become important. 



3.2 Correlation functions and relaxation time 

We investigate the temporal correlations in the system through a two-time single-site spin corre- 
lation function of the form 

n(4 -/- _l_ -A — £j=l s i(tw)Sj(t w + t) . . 

^ \?w ! T T,) — — jy ' [I J 

£*=1 «<(*!») 

Under equilibrium conditions this becomes a function of the relative time t only. The equilibrium 



correlation functions (averaged over 5 runs) for a range of temperatures are shown in Figure 7(a) 



As (3 increases, plateaux develop, revealing that again two-step relaxation is taking place. We also 



show examples of typical out-of-equilibrium correlation functions in Figure 7(b); one can see aging 
behaviour, with the correlation function showing dependence on the waiting time t w after a rapid 
quench from infinite temperature. 

Following the procedure in |p[,p7||, we may define a relaxation time r r as the time at which 
the equilibrium correlation function decays to e _1 . This is plotted on Figure ^ against inverse 
temperature; the data can be reasonably fitted by an Arrhenius curve of the form 

T r = Ae B ' T (8) 

where A, B are constants. The solid line superimposed on Figure || corresponds to A — 0.0166, B = 
2.535. This indicates that this model displays strong glassy behaviour, in agreement with the 
results from the topological model, where an offset Arrhenius law fitted the data considerably 
better than either a power law, or a Vogel-Fulcher law |Q. The anomalous upturn in the curve 
at very low j3 is due to the restriction that the spins may only take the values ±1 or 0; when the 
density of non-zero spins is very high, it becomes likely that some defects will sit next to each 
other, in locally 'stuck' configurations which slow the decay of the correlation function (see Figure 

However, if we look closely at Figure BL we see that at high inverse temperature the data 
seems to be drifting below the curve. If we consider more carefully the expected form of C(t) 
it becomes clear that the above definition of T r is not the most appropriate one, since again, as 
in the toplogical model, there are two decay processes. One can understand the origin of the 
plateaux by considering further the dominant processes involved in evolution of the system, which 
are analogous to those in the topological model |Q. In equilibrium, at low temperatures (high (3) 
there are few defects present in the system. The initial fast decay from C(t = 0) = 1 is due to 
dimers diffusing freely through the system, and thus moving the system away from the starting 
configuration. Isolated defects, however, need to either absorb or create a dimer in order to change 
position; this happens on a much longer time-scale than diffusion of the dimers. Therefore again we 
have two time-scales in the model - fast dynamics due to diffusion of the dimers, and slow dynamics 
due to movement of the isolated or stuck defects through absorption/creation of dimers. As noted 
earlier, the former is temperature-independent as dimer diffusion costs no energy; however, the 
position of the plateau and the subsequent departure from it is dependent on temperature. In 
particular, the depth of the drop to the plateau from the initial C(t = 0) = 1 is determined by 
the equilibrium concentration of spins of ±1 in local configurations that are free to move with no 
change to the energy (these local configurations include the dimers). One can therefore suggest 
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Time t (in units of N) 



(a) Equilibrium correlation functions C(t) for (from left to right) 
(3 = 1,2,3,3.5,4,4.5,5,5.5,6. 




Time t (in units of N) 



(b) Out of equilibrium correlation functions C(t w ,t w + t) for f3 = 6 for (from 
left to right) t w = 10 2 N, 10 3 N, 10 4 N. 



Figure 7: Correlation functions both in and out of equilibrium 
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nverse Temperature 



Figure 8: The logarithm of r r against inverse temperature. The solid line is an Arrhcnius 
law of the form r r = 0.0166 e 2E535 / T . 




Figure 9: 'Stuck' configurations. On the left, one sees that two spins of +1 adjacent to each 
other are unable to be flipped within the isolated cluster, regardless of what values the other two 
spins take (and similarly for two adjacent spins of —1). On the right, we see that two next-to- 
nearest neighbours of opposite spin cannot be flipped, regardless of the values of the other two 
spins. A favourable configuration external to the cluster can provide an opportunity to 'unstick' 
these cells. 
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Time t (in units of N) 

Figure 10: Correlation functions for, from left to right, (3 = 3,3.5,4,4.5,5,5.5,6. The solid lines 
superimposed are the best fits of the form C(t) = a e -t / Tl + (1 — a) e -t / T2 . 



that C{t) might be a sum of two functions in the following way: 

C(t)=af(t,n) + (l-a)g(t,r 2 ) (9) 

where a and r 2 are functions of T, but t\ is a constant, and a is the height at which one would 
expect to find a plateau. / and g are functions to be determined, but they must be monotonically 
decreasing functions of t, satisfying /(0, ti) = g(0,T 2 ) = 1 and /(oo,Ti) = g(oo,T 2 ) = 0. 

We tried to fit equation (||) to the data using exponentials for both / and g such that C{t) is 
of the following form: 

C(t)=ac- t / Tl + (l-a)e' t/T2 . (10) 

The results are shown in Figure [l^; in order to fit more accurately, the system size has been 
increased to N = 160000. The fit is extremely good at high values of f3 (low temperature), 
although at lower values there is some deviation. We also note that the fitted form drifts below 
the data at very low values of C (t) ; there may be some correction to this form which we have not 
taken into account. We do not claim that this form is exactly correct; nevertheless, it is a useful 
approximation that may allow us to separate out the two time-scales. Previously we stated that 



we expect t\ to be independent of temperature - Figure 11(a) shows that this seems to be the 
case at low /3, although as (3 increases it becomes harder and harder to fit T\ accurately due to 
the extremely high position of the plateau. We have shown the error bars on a few of the points 
to give some idea of the difficulty in accurately performing this fit at high (3; one can see that it is 
impossible to say anything sensible about the functional form of Ti{(3) for (3 > 4.5. However, the 
naive theory as discussed earlier in connection with E(t) gives Ti = 2 and the data is in general 
accord. 

We turn now to r 2 and consider the dominant processes involved in relaxation in the (3- 
relaxation regime. These are, as previously mentioned, absorption and creation of dimers. Creating 
a dimer costs 2 units of energy. Each dimer will rapidly diffuse freely through the system until it is 
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absorbed by a defect; this happens quickly as it is an energetically favourable process. Therefore 
we have energy barriers of 2 and probabilistic barriers of 2(3 in this regime. Absorption of pre- 
existing dimers also has a characteristic time-scale of ~ e 2 @ as the fraction of cells occupied by 
dimers scales as e _2/3 . This is reflected in the behaviour of the second time-scale as shown 



in Figure [11(b) : T2 exhibits the Arrhenius behaviour of T2 = Ae B @, with B = 2.12. We have 
shown error bars on a few of the points; these tend to suggest that the value of B is not exactly 
2; however, it is encouraging for the value to be so close given that this is a very crude theory. 

The plateau parameter a, is somewhat more complicated. As noted earlier, the initial fast 
decay of the correlation function is dominated by all the local configurations that can move freely, 
without any energy costs. Thus we expect the plateau height to occur roughly at a value of 1 
minus the fraction of total spins which can move freely in equilibrium conditions. This fraction 
contains the dimers, but also certain configurations of like-pairs (i.e. +1,+1 and —1, —1) which 



can oscillate, as shown in Figure 12, although they cannot delocalise without interacting with 



dimers. In equilibrium the probability that a given defect is part of a +1, —1 dimer is: 

M1M0) 2 = j^L- (ii) 

where p(s) is the probability of a spin of value s (the factor 6 comes from the fact the second defect 
may be situated on any of the 6 nearest neighbours to the original defect). The probability that 
any given defect is part of an oscillating like-pair is given by 12p(l)p(0) 2 (the factor 12 is because 
there are 12 positions at which the second defect may be situated to make up an oscillating like- 
pair). However, subsequent to t = one supposes that at any instant in time half of the oscillating 
like-pairs will be in exactly the same position as at time t = 0, and the other half will be in the 
alternative position. All other processes which cost no energy involve at least 3 defects, and thus 
are suppressed in comparison by factors of e^ 3 or e 2/3 . As we are investigating the region of (3 > 2.5, 
we neglect those, and simply suggest: 

Figure 11(c) shows this predicted curve of a against the simulation results for fitting with a sum 
of two exponentials, from which good agreement can be seen. Therefore for all three parameters 
(ti,T2,cx) the results of the simulations provide some support to the theory of the dominant 
processes involved in the evolution of this system. 

Note that the relaxation behaviour of the equilibrium correlation function C(t) differs from 
that of E(t) in two important respects. One is that for C(t) the system is always in macroscopic 
equilibrium so that the macroscopic distribution of non-zero spins fluctuates around a constant 
value throughout the dynamics, and all the observed results are due to re-arrangements of the 
location of these non-zero spins i.e. annihilations are balanced by creations on a macroscopic 
level. The second is that one has merely to move a non-zero spin in order to affect C (t) , whilst for 
E(t) to decay from the infinite temperature starting configuration one must actually annihilate 
±1 spins overall. 



3.3 Response and overlap functions 

We continue to show that this simple spin model behaves in the same fashion as the topological 
froth by studying response functions in relation to the fluctuation-dissipation ratio. Again we 
concentrate on the single-site (averaged) case, for which we require the linear response at a site 
to an infinitesimal perturbation field at the same site. The system is quenched from (3 — 1 to the 
temperature required, and then allowed to evolve at that temperature until time t w when a field 
of magnitude h and random sign Ci — ±1 is applied. Therefore the perturbation to the energy 
introduced by the field/charge interaction is: 

N 

AE(t)=h^e iSi (t)8(t-t w ) (13) 

i=l 
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Figure 11: The behaviour of the fitting parameters ti,x 2 and ex with temperature; 

obtained by fitting C(t) = a c~^ Tl + (1 - a) e'^ T2 . 
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Figure 12: Zero-energy moves for oscillating +1 like-pairs. An identical configuration exists 
for —1 like-pairs. 



where 6(t — t w ) is the Heaviside function: 9(t — t w ) = lVt > t w ; otherwise, h is a carefully 
chosen compromise which gives both linear response and a reasonable signal-to-noise ratio. The 
quantity which one measures is then the linear response function G(t w , t w +t): 

g(W»+«) = %5T%^ (14) 

One expects a parametric plot of — T G(t w , t w + 1) against C(t w , t w + 1) to have a slope of -1 where 
the equilibrium fluctuation-dissipation ratio is upheld. Breaking of this conventional equilibrium 
ratio is characteristic of aging in glasses f28| , |29[l ; the form of the slope when it is broken provides 
some information about the nature of the system. 

Figure |l3| shows the results at various temperatures: these show the same features as the topo- 
logical froth, namely a breakdown of the fluctuation-dissipation relation when t w is too short for 
equilibration to have occurred (for t w — > oo the system is already in equilibrium when the pertur- 
bation is introduced and FDT holds) and a reduction in the magnitude and an eventual change 
in sign of the slope as t increases and the correlation function decreases. The non-monotonicity 
is a consequence of the existence of an absorbing equilibrium state: competition exists between 
the field, which encourages the non-zero spins to settle on energetically favourable sites and thus 
increases the response, and the natural relaxation to equilibrium, which removes non-zero spins 
altogether thus reducing the response (one should recall that spins of value make no contribu- 
tion whatsoever to the response). At time t w (that is for t = 0) when the field is switched on, 
the response is on average zero. If t w is short enough that the fast processes remain dominant 
(i.e. E(t w ) has not yet reached the intermediate plateau) then these energetically favourable fast 
moves are quickly carried out. This simultaneously removes many non-zero spins, and also settles 
many non-zero spins on energetically favourable sites; the former process does not on average 
decrease the response since at t w it is zero on average anyway, and the latter process increases the 
response. Once the dimer concentration has relaxed to the equilibrium level (i.e. after the onset 
of the plateau in E{t w + t)) the response increases more slowly as the evolution of the system 
is dominated by the movement of isolated defects, which occurs on much longer time-scale. The 
applied field causes these isolated defects to tend to settle on energetically favourable sites; how- 
ever, as time passes they will be removed from the system until the equilibrium concentration is 
reached. Therefore one expects the response to peak before settling at a finite value once equilib- 
rium has been reached (the final value reached is dependent on the waiting time t w through the 
normalisation of both the response and the correlation functions). 

In order to distinguish between different types of aging, Barrat, Burioni and Mezard suggested 
the study of the overlap Q tu , (t) between two replicas . These replicas are identical configura- 
tions at time t w , but are subsequently evolved according to different stochastic thermal noise (but 
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Figure 13: Parametric plots of the response function against the correlation function 

(each averaged over 10 realisations of the charge distribution). 



of the same characteristic temperature), with Qt m (t) of the form: 



Qt w (*) = ^ 1=1 2n y ^ ( 15 ) 



where a is some order parameter of the system and the superscripts 1,2 refer to replicas 1,2. 
Barrat et al classify systems as either Type I or Type II models: for the former, the appropriate 
Qt w (t) (normalised to 1 at time t = 0) decays to a finite, non-zero value in the double limit 
lim^^oo lim^oo Q tm (t); this class includes models which are dominated by coarsening (for a 
review of coarsening see ]3l[]). For the class of Type II, Qt m {t) decays to zero in this limit; this 
class includes glassy systems M. In equilibrium (i.e. t w greater than the equilibration time of the 
system) there is no t w dependence so Qt„{t) = Q(t); one also finds that Q(t) = C(2t) (for details 
see @). 

For this model, the overlap we use is as given in equation ( ]l5| ) with = s^. We start from 
non-equilibrium conditions, where the system is quenched at t = from f3 = 1 to the temperature 
in question, and then allowed to run at that temperature until time t w when measurements 
commence. In Figure |l4|we show the overlap and C(t w , t w + 2t) against time for (3 — 5; initially Q 
is almost identical to C(t w , t w + 2t), but it drops below C(t w ,t w + 2t) at longer times. For larger 
values of t w this takes longer to happen since these systems start off closer to equilibrium. We 
also show the equilibrium curves, for which Q(t) — C(2t). 

Given that Q(t) = C(2t) in equilibrium, we can check our proposed functional form for the 
equilibrium correlation function through a parametric plot of the overlap against C'(t). If C(t) is 
indeed of the form given in equation (|10|) then one expects: 

Q(t) = C{2t) = a e- 2t/Tl + (1 - a) e~ 2t/T2 (16) 

Thus for long times one would expect to find: 

Q(t) ~ (1 - a) e-»/« ~ (17) 

(1-a) 
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Fi gure 14: The behaviour of C(t w ,t w -\- 2t), Qt w i^) with time for j3 — 5. Qt^(^) is given by 
the dashed curves and C(t w ,t w + 2t) by the solid curves. From lower curve pair to upper curve 
pair, t w = 107V, 10 2 7V, 10 3 iV, 10 4 7V. The uppermost pair is that of equilibrium, with Q(t) = C(2t). 



and for short times: 



Q(t) ~ a e- 2 */- 1 + (1 - a) ~ { ° {t) + - 1} " + (1 - a) (18) 



n 



Figure jig shows parametric plots of the overlap Q(t) against the correlation function C(t) (in 
equilibrium) for different values of f3. These plots should be 'read' from the top right corner i.e. 
t = occurs when Q(t) — C(t) — 1, and long times correspond to low values of Q(t),C(t). We 
have plotted on each the expected short time and long time behaviour as given in equations ( |l7j ) 
and (|8|), where the values of a are those fitted in the previous section. One sees that the short 
time expression fits the data extremely well, supporting our hypothesis that the initial decay of 
the correlator, from 1 down to the plateau, is exponential. The longer time behaviour initially fits 
very well for both (3 — 4 and 5 (Figures 15(a)| , |l5(b)| ), although as time goes on the theoretical 



curve drops below the data. This tends to suggest that at long timc-scalcs there is some correction 
to this fitted form which we have not taken into account; the second relaxation may be some kind 



of modified exponential rather than pure exponential. For the case of (3 = 3 (Figure 15(c)) the 
values of Ti,T2 are much closer together and thus the two time-scales are not so well separated. 
Therefore one does not see a well-defined plateau in the equilibrium correlation function (see 



Figure 7(a)), and the cusp in this parametric plot is also not as clear. The cusp in these plots 
is a direct result of the existence of the plateau; since Q(t) = C(2t) in equilibrium, the overlap 
function Q{t) reaches the plateau before the correlator C(t). Thus there is a time period for which 
Q(t) is effectively stationary whilst C{t) is still dropping fast. This is followed by a time period for 
which the plateaux in both functions overlap, and therefore both are stationary, and then there 
is a regime in which Q(t) drops away from the plateau whilst C(t) is still stationary. This tells 
us that whilst the dominant process is diffusion of the dimers, the two copies of the system are 
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Figure 15: The overlap function Q(t) against the correlation function C(t). In each case 
the dot-dashed curve is the expected short time behaviour and the dashed curve is the expected 
long time behaviour, if equation (FRl) holds. The values of a are those fitted in the previous section. 



restricted to a narrow area of phase space; this is because the isolated defects have not yet moved 
in either copy, and thus the overlap will be high. It is only once the activated processes become 
dominant that the two copies can move well apart from each other. 

Figure |l6| shows a parametric plot of the overlap against the correlator for the non-equilibrium 
case i.e. Qt w (t) against C(t w , t w + t). The initial behaviour is independent of t w and in fact follows 
the short time behaviour we expect in equilibrium. However, as C{t w ,t w + 1) drops below 0.9, the 
curves fall below the equilibrium behaviour and we see evidence of dependence on t w . The cusp 
that develops is sharper for larger t w , with the behaviour tending towards that of equilibrium. For 
the smaller values of t w there is no channelling in phase space because the starting configuration 
is well away from equilibrium, and thus there are many different energetically favourable routes to 
be taken. It is clear from these figures that our model falls into the class of Type II, as the overlap 
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Figure 16: Qt w (t) against C(t w , t w + t) for /3 = 4. From lower to upper curve, t w = 
ION, 10 2 iV, W 3 N, 10 4 7V. The dotted curve is that of expected short time behaviour for Q(t) 
against C{t) in equilibrium. 



decays to zero rather than a finite value as it would do for Type I systems; this is as expected as 
we believe our model to be glassy rather than dominated by coarsening. 

To summarise our findings so far, we have shown that the D > model docs indeed give results 
that are quantitatively similar to those of the topological model, whilst having the advantage of 
being simpler, computationally faster and more suited to analytic study. We have developed a 
conceptual picture involving both fast and slow dynamics, and many of the features of this model 
can be described in terms of this picture. We find good agreement between theoretical predictions 
and data for the behaviour of the correlation function, energy and overlap function. 

4 D < 

4.1 Relaxation dynamics and correlation functions 

We shall now turn our attention to the model with D < 0, setting D = — 1. Unlike the previous 
case, this model does not have a unique absorbing ground state; instead, there are a great many 
degenerate ground states, although some of these are not accessible using our dynamical rules. This 
raises the possibility that in preparing an 'equilibrium' system by randomly placing throughout 
the lattice the correct number of ±l's for that temperature, one might pick out an inaccessible 
configuration. However, the probability of this occurring is so small as to be negligible, and we 
have checked that the results obtained by this method do not differ from those obtained through 
waiting long enough for equilibration to occur. Thus when we refer to an equilibrated system, we 
mean one that has been prepared through a random allocation of the correct number of non-zero 
spins. 

We expect the D < model to also show glassy behaviour with two different time-scales, but 
as the system now favours s, = ±1 rather than Sj = 0, the equivalent of the free-moving dimers 
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(a) The behaviour of the energy under cooling (D < 0). The values of t w are the 
waiting times at each point. 




(b) The behaviour of the energy after a rapid quench (D < 0). The values of t 
are the times, subsequent to the quench, at which the energy is measured. 



;ure 17: Energy against temperature for slow cooling and rapid quench. 
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(a) Energy against time. 
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(b) Energy against T lnt, where t is measured in units of N. 

Figure 18: The behaviour of the energy with time. 
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(a) The persistence of non-isolated zero-spins for (b) The persistence of non-zero spins with at least 

D < 0. one opposite neighbour for D > 0. 

Figure 19: Persistence functions for both D < and D > 0. 

will be fast-moving pairs of spin zero, and the analogue of the slow-moving energetically-trapped 
defects will be isolated zero-spins. With the correct choice of observables, one expects to see all 
the same features as with the D > model. However, although one expects similar qualitative 
behaviour,^ the quantitative behaviour should be different - this is because A + A — ► processes 
replace those of type A + B — > and also because the zero-spin dimers cannot move quite so easily 
through the ±1 background. As shown in Figure ^| there are certain configurations that simply 
cannot move. 



Simulations were again performed for N = 9900. Figure |l7(a)| confirms that we do see glassy 
behaviour when the system is cooled at different cooling rates. A plot of the energies attained 
after running for a variety of times t w at various temperatures T from a starting configuration 



corresponding to infinite temperature are shown in Figure 17(b) ; the presence of activated pro- 



cesses is indicated by the clear minima and by the plateau a,tjf^ —0.94 below which the system 
cannot penetrate at low temperatures, even after the longest waiting times. Again one can see this 
plateau clearly in a plot of energy against time subsequent to a quench from infinite temperature 



appear (Figure |l8(b)|) 



(see Figure |18(a)D; on e can again re-scale the time axis to T In t in order to see the staircase shape 



Because the background does not necessarily allow the dimers to move freely, one must consider 
the possibility that the dimer diffusion may be dependent upon the density of non-zero spins. In 
equilibrium, this density is dependent upon the temperature; therefore we have investigated the 
persistence of non-isolated zero-spins under equilibrium conditions. In order to do this, one can 
identify all zero-spins that have at least one neighbouring spin which is also zero in the starting 
equilibrium configuration; as the system evolves, one can measure the fraction of these that have 
NOT been involved in a move. Figure 19(a)| shows these results for D < 0; one can see that 



the persistence is weakly temperature-dependent. One cannot sensibly investigate the effect at 
lower temperatures because there are so few zero-spins present in equilibrium conditions. For 
comparison, we also show the results for the D > model, where the persistence is defined as 
the fraction of non-zero spins with at least one neighbour of the opposite sign that have NOT 



undergone a move. The results are shown in Figure 19(b) and show no temperature-dependence 



This result has implications for the form of both the correlation function and the energy. 

^except in pathological cases 
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(a) The energy for (from left to right) (3 = 5, 6, 7 fitted with equation (|19|); the 
fits are the dashed curves. The parameters (ai, mi, K\, K2) are as as follows: for 
P = 5, (-0.942,1.02,0.49,0.52); for f3 = 6, (-0.941,0.93,0.50,0.51) and for (3 = 7, 
(-0.941,0.93,0.49,0.54). 



-0.65 




Time t (in units of N) 



(b) The energy fitted with equation (pcj); the fits 
are the solid lines. (01,02, mil Ki, rri2, v, K2) are as fol- 
lows : for 13=5, (-0.936,0.61,0.748,0.98,0.0385,0.638,0.62); for /3=6, 
(-0.937,0.57,0.77,0.977,0.0389,0.64,0.59) and for (3=7, 
(-0.937, 0.66, 0.46, 1.0, 0.0313,0.474, 0.59). 



Figure 20: Fits to the energy. 
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In order to produce a fit for the energy relaxation from a starting configuration corresponding 
to infinite temperature, one must think very clearly about the processes involved. The initial fast 
decay to the plateau involves two pairs of zero-spin dimers annihilating; this is of type A + A — > 0, 
which one expects to give a("5 dependence in the energy in the asymptotic limit |L8],|li| 0. At 
this stage the density-dependence of the dimer diffusion is not likely to be strong enough to cause 
substantial deviation from the A + A — > behaviour so one expects the initial decay to behave as 
(1 + mii) ; however, it is not trivial in this case to determine what the value of mi might be. 
Because the zero-spin dimers are not completely free to move diffusively through the background 
we are unable to produce an accurate estimate; we simply expect a value of vn\ of the order of 1. 
For the slow decay from the plateau to equilibrium, the density-dependence of the dimer diffusion 
will have a substantial effect because there are far more ±1 spins present as the system flows closer 
to equilibrium. The slow process involves the pairing of isolated zero-spins through the mechanism 
of dimer creation and absorption; it is unlikely to behave exactly as an A + A — > process given 
that the dimers cannot diffuse freely to facilitate it. One can fit this latter part separately, and 
one finds it takes the form a\{l + ie -2 ^)~ K ; the factor of two in the exponential was fitted as a free 
parameter and a result of almost exactly 2 was obtained for every temperature. This is in keeping 
with the energy barrier of 2 involved in creating a dimer. The parameter a\ is naturally associated 
with the plateau value and clearly takes a value ~ —0.94. The value of k is approximately 0.6; this 
is substantially slower than the behaviour one would find asymptotically with a pure A + A — > 
process and is due to the inhibited movement of the dimers. 

Having fitted the latter part, we attempted a fit of the full dataset of the form: 



E 
N 



f-~ -Oi J (1 + mit) Kl + (oi - e eq ) (1 + te' 20 ) K2 +e eq (19) 

where a%, Ki, K^, mi are all parameters to be determined, and e eq is the equilibrium energy per cell 
at the temperature in question. One cannot fit this well to the data, as shown in Figure I'm.'! 



the decay is faster than that of a power law in the latter stages of the decay to the intermediate 
plateau. Thus we must think again about the processes involved in the relaxation of the energy 
to the plateau. 

As mentioned earlier, besides the A + A — > fast processes, there is also a fast process involving 
a dimer interacting with a defect to leave an isolated defect - this is of type A + C — > + C, and 
thus typically gives a stretched exponential for the asymptotic behaviour of the energy with time. 
We did not need to include these in the fit to the energy for D > 0; however, it is clear from the 



poor fit in Figure 20(a) that in the D < case we cannot neglect them. We therefore expect the 



energy density E/N to be approximated by the following form: 

f = ("5 " °i) ("2 (! + m ^ Kl + C 1 - «2)e- (m2t) ") + (ai - e eq ) (l + te^)"" 2 + e eq (20) 

where a\, K2 can be fitted separately from the decay from the plateau to equilibrium, and 02, Ki, 
v,m\,m2 are parameters to be determined. The theory would suggest that Ki should be close 
to 1, and that mi should be of the order of 1. Figure |20(b)| shows these fits superimposed on 
the data; it will be noted that agreement is excellent. The values of the parameters are given in 
the caption; in particular, one notes that K\ is extremely close to 1 in each case, and that mi is 
indeed of the order of 1. We have already discussed the power law decay from the intermediate 
plateau to equilibrium. Without developing a much more complex theory one can say little about 
the parameters in the stretched exponential term. 

With regards to the correlation function, in order to study the D > model we chose a function 
which focussed on the defects. We do the same in this case, and thus the correlation function we 
measure is: 

C(t w ,t w+ t) = - Peq (0)) (21) 

t-Peq{U) \ Ei=lVt ro ),0 / 



4 In fact the critical dimension d c = 2 for this theory, and one expects logarithmic corrections at this point, but 
we shall ignore these. 
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Figure 21: Correlation functions in equilibrium conditions. From left to right, (3 - 
2.5,3,3.5,4,4.5,5. The superimposed fits are of the form: C(t) = ae^ 1 + (1 - a) e"^/^)" 



where p e q{0) is the equilibrium density of zero spins at the temperature of interest. The somewhat 
unusual normalisation is necessary to allow the correlation function to decay to zero rather than 
to a finite value equal to p eq (0); this finite plateau comes about because a fraction p eq (0) of the 
cells that have spin zero at time t = can be expected to have spin zero at any later time. In the 
D > model, we did not have to take into account this effect since the contribution from spins 
that remain +1 (or -1), and those that swap from +1 to -1 (or -1 to +1) cancel out. Figure pi] 
shows that this function does indeed produce similar results to those of the D > model: again, 
we can see clear evidence of two-step relaxation. The fits superimposed upon the data are of the 
following form: 



C(t) = a c-*/ Tl +(l-o) e^/^'' 



(22) 



It will be noted that these fits are extremely good, and also that the relaxation in the long-time 
region follows a stretched exponential rather than an exponential as in the D > case. This is 
a result of the temperature dependence of the dimer diffusion, which slows the decay from the 
plateau. Thus 7 in Figure |22(d)| is less than one in all cases. We do not have a full theory of the 
behaviour of the system in this region and thus cannot comment further on the behaviour of this 
exponent. 

One might naively expect the predicted value of t\ to be altered by the fact that the zero-spin 
dimers are not necessarily free to move through the ±1 background. In fact this is not the case: 
Ti is the time-scale for the dimers that are free to move only, and we can still expect those free 
dimers to move with a time-scale of 2, independent of temperature, exactly as in the D > 
case. Figure [22(a)| shows the values of n against inverse temperature obtained from fitting the 
correlation functions with equation (^2|) - this data is in keeping with a temperature-independent 
value of ri = 2. Those dimers that are not free to move do not contribute to the initial fast 
decay of the correlation function. The parameter that the extra jamming does alter is a: one 
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Figure 22: The behaviour of the fitting parameters ol,t^,t-2 and 7 with temperature. 

The data was obtained by fitting the form C(t) = a e _t / Tl + (1 — a) e^*/ 1 " 2 ' 7 to the equilibrium 
correlation functions. 



expects the correlation function to decay to a plateau value which is one minus the density of free 
zero-spins (although one must remember to normalise correctly as in equation Thus we 

have to calculate the probability of obtaining a zero-spin dimer which can move (shown in Figure 
4(c) , but one should now think of the ±l's as being the background and the zero-spins as being 
the dimer), and also of obtaining a zero-spin pair which can oscillate (as in Figure 12). This gives 
a probability of 24p(0)p(l) 2 ; thus after normalisation, we expect a to behave as: 



12e" 



-0 



(1 + 26^)2 



(23) 



Figure 22(c) shows this curve superimposed upon the fitted values of a against temperature; the 
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Figure 23: Parametric plots of the response against the correlation function for differ- 
ent temperatures and waiting times. 



data shows reasonable agreement with the theory. 

One can see from Figure 22(b) that we again observe Arrhenius behaviour for T2. In this 
case, however, the best-fit Arrhenius law is ti ~ e 2 ' 4/3 , whereas the energy barrier argument would 



suggest T2 ~ e ' . As mentioned earlier, we lack a full understanding of the behaviour in this region, 
and can only say that this discrepancy is probably also due to the fact that the ±1 background 
does not allow the zero-spin dimers to diffuse freely. 



4.2 Response and overlap functions 

We continue our study of the D < model by observing the response of the system to an applied 
field. As before, a field hti is applied at time t w , with e, = ±1 randomly at each site. The 
perturbation this introduces to the energy is AE(t) = €i$s(t),o ^(^ — ^w)i where 9(t — t w ) 

is again the Heaviside function. The observable we then measure is the following linear response 
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(a) (3 = 2.5 (b) (3 = 3.5 

Figure 24: Overlap Q(t) against C(t) in equilibrium. In each case the dashed line is the 
expected short-time behaviour and the dot-dashed line is the expected long-time behaviour, if 
equation (p2|) holds. The values of a, 7 are those fitted in the previous section. 



function G(t w , t w +t): 



Again one expects a parametric plot of —TG{t Wl t w + t) against the correlation function (defined 
as in equation (|2~l])) to yield a slope of -1 where the conventional FDT is upheld. 

Figure ^ shows such a parametric plot for a variety of temperatures and waiting times t w . 
One can see the main features of the equivalent D > plots in these figures: FDT is upheld for 
a time which increases as t w increases, and sometime after it is broken the curves display non- 
monotonic behaviour. As in the D < case, this turnover in the response after long times (i.e. 
low values of the response and of C(t w , t w +t)) is due to the isolated defects being eliminated from 
the system, thus decreasing the response. Note that although this turnover can been seen clearly 



in Figure 23(a) , one cannot always run the simulations for long enough to observe this effect at 
low temperatures and large waiting times; in fact, it is not observable for any waiting times for 
(3 = 6 (Figure |23(c)| ). However, in addition to this behaviour, at very low temperatures and for 



short t w we see an intermediate 'hump' appearing before the turnover due to the isolated defects 
see Figures E3(b)| , |23(c)| . This is a consequence of the presence of temperature-dependent dimcr 



diffusion; blocked dimers exist which take some time to become mobile and cannot diffuse freely 
through the system. Instead of being eliminated before they can make a substantial contribution 
to the response, some dimers persist and considerably increase the response before they are finally 
removed. This accounts for the intermediate humps shown for short waiting times in Figures 
|23(b)| , |23(c") ; at longer waiting times this effect is imperceptible because the dimers have moved 



closer to equilibrium. It is also imperceptible at higher temperatures (lower 0) because the dimers 
can move more freely and equilibrate more quickly. 

Let us turn now to the overlap function, defined in this case as: 

E^x^.oW 

Recalling that in equilibrium one finds Q(t) = C(2t), we can use the equilibrium overlap to test 
of our proposed form of C(t), as given in equation (|22|). If this equation holds we expect to find 
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that for short times: 



Q(t) ~ (ae- 2 */- + (1 - a)) ~ + (1 _ (26) 

and for long times: 

Q(t)~(l-a)e<-*/^~-^L. (27) 

(1 - a) 

Figure [24| shows the equilibrium results for f3 — 2.5 and 3.5; the superimposed curves are the 
expected short and long time behaviour using the values of a, 7 obtained from fitting the correlation 
functions with equation (|2^). The theoretical behaviour clearly fits the data very well, lending 
further support to equation (|2^) as a description of the behaviour of the equilibrium correlation 
functions. We also see that Q(t) decays to zero, thus placing the D < model in the Type II class 
along with the D > model. 



5 Concluding remarks 

We have studied a simple lattice-based spin model which has a non-interacting Hamiltonian, 
but constrained dynamics, and find it to exhibit both glassy behaviour and behaviour typical of 
diffusion-limited reaction models. A single parameter D distinguishes two types of ground state. 
By choosing D > 0, one can study a system with a unique ground state, evolving by way of a 
number of annihilation-diffusion processes which arc either fast temperature-independent or slow 
temperature-dependent diffusive processes, with the latter slower by a factor exponential in inverse 
temperature. We can categorise the fast processes as: A + A — > 0, A + A (or A + A) — > C + C, A 
(or A) + C (or C) — > + C (or C); and the slow processes as: C + C — > A (or A), where A and 
A are dimers and anti-dimers, and C and C are isolated defects of opposite sign. The isolated 
defects move isotropically, but the dimers and anti-dimers move anisotropically, and come in three 
different 'flavours' according to their orientation. The different flavours can also scatter amongst 
themselves via the process A a + A a — > A" + Al 3 , where a, (3 label different flavours. Since A ^ A 
and C ^ C, A + A processes are equivalent to the usual A + B processes, and C + C are equivalent 
to C + D processes. In this paper we study the full set of processes simulationally but only provide 
a simplified adiabatic theoretical fit. For D < 0, the ground state is highly degenerate and the 
system evolves according to the fast annihilation-diffusion processes ^4 + ^4^0, A + A — ► C + C, 
A + C — > + C, and the slow diffusive process C + C — ► A, where A corresponds to a pair of 
zero-spins and C to an isolated zero-spin. Again there are 3 different flavours of A corresponding 
to the three different orientations, and these can scatter through A a +A a — > A^ + A^. In this case 
the movement of A' a are hindered by the background due to the ground state degeneracy, and 
even the fast processes have some temperature dependence. As before we have studied the full 
dynamics simulationally but only provide a simplified adiabatic analysis. There is clearly scope 
for providing a full analytic theory. 

For both the D > and the D < case one finds two-step relaxation, on two different time- 
scales which are separable and can be attributed directly to the different processes. We find that 
the energy in the D > case can be fitted with the sum of two terms, each behaving like the 
asymptotic predictions of an A + B — > theory (for the slow processes, C + C — ► A, but the A' a 
are eliminated on a time-scale which is negligible compared to that of the slow process, so this 
behaves like C + C — > 0) . In the D < case one cannot fit the energy adequately in the approach 
to the intermediate plateau without including induced dimer absorption (A + C — > + C) along 
with the A + A — > fast diffusive processes. We also find that the slow diffusive process does not 
behave like a pure C + C — > process; this is in part due to the fact that the dimer diffusion is now 
temperature dependent. We have studied the correlation functions for D positive and negative; 
in both cases, a naive theory gives a predicted form for these which fits the data extremely 
well. Studies of the overlap in equilibrium serve to reinforce these results. An investigation of 
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the response function in both cases yields non-monotonic response curves, and for D < the 
temperature-dependent diffusion leads to more complex results: one finds intermediate humps for 
very low temperatures and short waiting times, which can be understood within the framework 
of the processes we have already discussed. Non-monotonicity has also been observed in the 
response functions of many other models which involve activated processes |p[-p|,p^| . To distinguish 
between Type I (coarsening) and Type II (glassy) tendencies we have examined an overlap function 
measuring the temporal auto-correlation of two independently evolving clones of a configuration. 
This demonstrates that the present system is of Type II, for D both positive and negative. 

In this paper we have employed a hexagonal basis for the cell edges. This is naturally motivated 
by analogy with a two-dimensional froth. It also corresponds to the case of the simplest non-trivial 
vertices, which have valence three, and consequently is special in that any cell has two nearest- 
neighbour cells which are nearest-neighbours of one another. Extensions are clearly possible, both 
to higher valence vertices in two dimensions and to minimal and non-minimal vertices in higher 
dimensions, but we do not pursue them here. We merely note that systems with valence greater 
than (d + 1), where d is the dimensionality, are more prone to sticking. 
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